Model for the resistance force acting on circular bodies in the imminence of rolling 
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The laws of friction are reasonably well understood for the case of blocks in contact with rough 
plane surfaces. However, as far as bodies with circular sections are concerned, the physics of friction 
becomes more involving and it is not possible to adopt a simple conceptual framework to explain all 
phenomena. In particular, to the best of our knowledge, there is no approach so far to the problem 
of the resistance force that opposes to circular bodies that remain at rest while acted upon by small 
external forces. Here we fill this gap by introducing a purely mechanical model based on both the 
elasticity theory and Hertz contact mechanics. Our approach furnishes a quantitative expression for 
the critical force beyond which rest can no longer be maintained. Besides confirming the expected 
proportionality of the resistance force with the load, our result contains no free parameters and is 
expressed solely in terms of physical properties of the problem, such as the pressure of the body per 
unit of superficial area, a relation between the Young modulus of the surface and its Poisson ratio, 
and the symmetry of the contact. 

PACS numbers: 62.20.Qp,81.40.Pq,46.55.+d,45.20.-d 
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I. INTRODUCTION 

When two solid bodies in contact move relatively to 
one another (or when a resting body is acted upon by a 
tiny applied force) there appears a resistance force, con- 
trary to the (onset of) motion, which we generically call 
friction. Estimates point out that up to 1.4% of the GNP 
of a country may be wasted with energy losses associated 
with friction [l( . Such a contact interaction permeates a 
number of aspects of the human life, being important 
even for everyday acts such as walking and wearing. De- 
spite the undisputed relevance of friction and its longterm 
history in science, it is fair to say that its basic physics 
is not yet fully understood. 

The first attempts to describe the physics of friction 
date back to the efforts of da Vinci, who empirically de- 
rived the basic laws of friction (posteriorly rediscovered 
by Amontons), namely, that friction is proportional to 
the load and independent of the area of the sliding sur- 
faces. Later on, many others contributed to the foun- 
dations of the science of Tribology. Eulcr proposed that 
friction was caused by interlocking irregularities between 
the surfaces. In addition, he introduced the coefficient of 
friction /i and was the first to make a distinction between 
static and kinetic friction. It is due to Coulomb the laws 
according to which there exists a maximum static fric- 
tion proportional to the load and the kinetic friction is 
independent of the sliding velocity. Later on, Bowdcn 
and Tabor [2j gave a physical explanation to the afore- 
mentioned laws of friction. According to their approach, 
(i) the contact between surfaces is due to asperities and 
(ii) the actual contact area, which friction ultimately de- 
pends on, is small compared to the apparent area. Thus, 
the larger the applied load the larger the real contact 
area, which in turn implies larger friction. 

As far as bodies with circular sections (mainly cylin- 
ders and spheres) are concerned, the physics of friction 



becomes more subtle. It has been recognized since long 
ago that in situations involving the so called rolling fric- 
tion it is no longer possible to neglect the deformations of 
the bodies [3J]. In fact, studies on this subject often ap- 
pear attached with the Contact Mechanics Theory [J-Q , 
whose founder father is Hertz. Although his formulation 
neglects adhesion and Van der Waals interactions, the 
contact area can often be analytically derived in terms of 
the physical properties of the bodies. The adhesive inter- 
action and the Van der Waals interaction are sequentially 
added to Hertz's approach via the JKR and DMT theo- 
ries, respectively (7|- Tabor demonstrated that both the- 
ories are successful, each one in a specific regime, which 
is defined by the materials and geometry involved. The 
gap between those two is filled by the Maugis-Dugdale 
theory. Also, it is worth noting that the JKR and DMT 
theories attribute the energy loss to hysteresis, not slip, 
what explains why it is barely affected by lubricants [8|. 
The reader is referred to Refs. 0, [l(| for a survey of all 
the aforementioned theories and their regimes of validity. 
In many treatments, rolling friction appears as a torque 
pP opposed to the rotation of the body, where P is the 
load and p is the coefficient of rolling friction. The lit- 
erature points out that p is in the range 10 -2 to 10~ 3 
11 1, although some discordance exists [12 1. 



The problem of rolling friction was investigated in a 
variety of interesting situations, most of them involving 
macroscopic regimes |13H2ll |. Recently, however, new in- 
sights have been put forward by adopting a microscopic 
perspective (see Ref. [22j for a review on this topic), by 
means of which the origins of friction arc explored in 
the atomic domain. Efforts in this direction, often con- 
jugate to friction- and atomic- force microscopy [ 23 M 25|. 
define the area referred to as Nanotribology J26H28J |. In 
this context, the problem of friction between two sub- 
strates separated by a thin lubricant film was theoret- 
ically addressed via molecular dynamics simulations in 



simple two-dimensional models involving Lennard- Jones 
potentials [29(. Experimentally, setups based on atomic 
force microscopy have been used (i) to directly measure 
the adhesion and rolling-friction forces between silica mi- 
crospheres [3(j] and (ii) to measure and control the rolling 
and sliding of carbon nanotubes at nanometric scale [3l| . 
This paper aims to contribute to an issue which, to 
the best of our knowledge, remains open so far. It refers 
to a very basic problem with rather nontrivial solution. 
We are concerned with the resistance force that acts on 
a circular macroscopic body (sphere or cylinder) while 
the static equilibrium is maintained, i.e., while the body 
remains at rest. This is analogous to the situation in 
which a block remains at rest on a rough horizontal sur- 
face while acted upon by a tiny force F = Fe x ; a fric- 
tional force f = —Fe x appears in the contact area which 
prevents the motion to occur. This certainly cannot be 
the whole story for circular bodies. Indeed, it is easy 
to see that if the external force is applied horizontally, 
e.g., towards the center of mass of the body, those forces 
would produce a resulting torque in relation to center of 
mass and the static equilibrium would no longer exist. 
Nevertheless, we do know from our everyday experience 
that static equilibrium must exist for a sufficiently small 
F [33. Our approach consists of using results of the 
Hertz Contact Mechanics to build, out of a minimum set 
of assumptions, a mechanical model for the resistance 
force acting on macroscopic circular bodies. We derive 
an upper bound for the maximum force that can be ap- 
plied before the body starts to move. Free from any 
fitting parameter, our formula is shown to depend only 
on physical quantities, such as the load, the geometry of 
the body, and the elastic properties of the surface. 



II. MODEL 

Our aim is to devise a mechanical model without any 
appeal for microscopic details of matter, such as its 
atomic structure. As mentioned above, we are concerned 
with the forces that act while the static equilibrium is 
preserved. Moreover, we focus on cases in which the 
static regime is necessarily replaced by rolling without 
sliding. For concreteness, we can imagine the onset of 
macroscopic motion of a bowling ball (initially at rest on 
a horizontal surface) when acted upon by a force that 
smoothly increases from zero. 



f = — fe x (friction) and N 
and N = P. 



Ne y (normal), with / 



F 




FIG. 1. (a) An external force F = Fe x , smaller than the 
maximum static friction, is applied on a cubic block of edge 
21 and weight P = ~Pe y initially at rest. Contact forces 
f = —fe x (friction) and N = Ne y (normal) maintain the cube 
at rest. The rotational static equilibrium is preserved when 
x = §1(1 + tan#). (b) The two contact forces are merged 
into a single contact force C, whose inclination is given by 
tana = P/F. 



Now, in order to preserve the rotational static equi- 
librium, the normal force has to act a distance x = 
§1(1 + tan#) far from the center of symmetry [see 
Fig. [Ha)], a result that immediately follows from the 
basic laws of classical mechanics upon the assumption 
that the total torque vanishes. We see that the point 
of application of the normal force moves sideways as ei- 
ther F or 9 increases. In practice, this effect may de- 
rive from unperceivable rotations and deformations that 
make the microscopic contact area move laterally. When 
F(I + tan 6) > P, one has that x > I, implying that 
the rotational equilibrium can no longer be maintained. 
The point of application of the friction force f does not 
influence the torque in relation to center of mass, but pre- 
sumably it should be where the contact is more effective. 
One may naturally expect that both f and N effectively 
act on the same region. 

This observation makes it possible to merge the friction 
and the normal forces into a single contact force C [see 
Fig. HUb)], from which we may derive f = (e x ■ C)e x 
and N = (e y ■ C) e y . Accordingly, the static equilibrium 
will occur for C + F + P = and x = §1(1 + tan0). Of 
course, we are preserving the idealization of a single point 
of application for the contact force, which actually is a 
resulting force that encapsulates the global reaction of 
the surface. This approach is particularly interesting for 
problems involving circular bodies, for which the effective 
contact area is rather smaller. 



A. Contact force 

Let us start by considering the problem of a cubic block 
of weight P = —Pe y and edge 21 acted upon by a force 
F = Fe x whose modulus is smaller than the maximum 
static friction fj,N. The situation is depicted in Fig. [T] 
As usual, to maintain the translational static equilibrium 
one should demand the contact forces to be given by 



B. Rigid bodies and elastic surfaces 

We focus on regimes in which an initially plane surface 
deforms elastically under the weight P of an ideally rigid 
circular body of radius R. The situation is illustrated in 
Fig. [SJ When the body is acted upon by a sufficiently 
small external force F, a contact force C appears in the 



point of application O in order to preserve the static equi- 
librium. The angular displacement of O from the vertical 
is denoted by ^ and the penetration length by d. From 
the geometry of Fig. [2] it is straightforward to check that 
# e [0, * g ], where 



* S = 



1/2 



(1) 



"g" standing for "geometric." The quality of this ap- 
proximation is determined by the elastic properties of 
the surface, in particular by its Young's modulus, E. The 
larger the value of E the smaller the deformation d and 
the better the approximation ([T]). Note that the very 
notion of ty g is not much accurate because the actual 
contact area is smaller than that suggest by the geom- 
etry. In other words, because of the smoothness of the 
surface — there should not be any discontinuity in the first 
derivative at any point of the surface — the last point of 
contact actually occurs at an angle slightly smaller than 
the geometrical value given by Eq. (JXJ) . 




FIG. 2. A rigid circular body of radius R and weight P de- 
forms an elastic plane up to a maximum penetration length 
d. Under the application of a tiny force F a reaction force C 
emerges which guarantees the static equilibrium. The angular 
position of the point O at which C applies is *?. 

An important advantage of the present framework is 
that the Contact Mechanics Theory furnishes several an- 
alytical results describing the characteristics of the elas- 
tic deformation in terms of the physical properties of the 
bodies involved [J, Q . A fundamental building block of 
this theory is the formulation of the correct pressure dis- 
tribution. Theoretical treatments of these problems can 
be found in Refs. [J, la, 0, l33| . Of particular interest for 
our purposes are the formulas of the penetration length 
for a sphere of radius R, 
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phere 
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and a cylinder of radius R and length L, 

4 P 



"cylinder 
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(2a) 



(2b) 



As usual, the results are written in terms of the parame- 
ter k = E/(l — v 2 ), here referred to as medium constant, 
as it gives a relation between the Young modulus E and 



the Poisson ratio v of the medium. Although the re- 
sult (f2"a|) is universally accepted, for the cylinder there is 
still some discordance. Indeed, by direct inspection one 

may verify that 7 = § In (^p- 5 ^-) m Ref. 33( whereas 
7 = 1 in Re f. |4j , the former enjoying some experimen- 
tal support [34| . In addition, it is worth noting that 
these results accurately hold only in the Hertz regime, 
i.e., when the adhesion forces can be neglected in com- 
parison with the weight of the body. Off this regime, 
the above formulas should be superseded via appropriate 
theories [lfj. Our model, however, will keep applicable, 
provided we include the adhesion in the balance of forces 
and torques. 



C. Maximum resistance force 

For blocks, while the static equilibrium is maintained 
under the action of an external force, the horizontal re- 
action of the surface is generically called static friction. 
For circular bodies there is some ambiguity, for some au- 
thors associate the term with the situation in which the 
body rolls without sliding. To avoid confusion, here we 
call resistance force the total horizontal reaction offered 
by the contact force. Formally, however, we preserve the 
notation 



f = (e x ■ C)e x , 



(3) 



which obviously leads to f = —Fe x in the present prob- 
lem. It is immediate that the resistance force is null for 
F = 0, as in this case C = Pe y . If we fix the modulus F 
of the external force and increase the angle of application 
9, then \1/ should increase to balance the torque, while the 
modulus of C keeps fixed to guarantee the translational 
static equilibrium (see Fig. [5]). On the other hand, if we 
fix 9 and increase F, then both the modulus and the in- 
clination a (defined as in Fig. [1} of the contact force in 
relation to the x axis has to vary in order to ensure that 
C + F + P = 0. In this case as well, "J has to increase 
to preserve static equilibrium. Clearly, VP is a function 
of the applied torque, i.e., ^ = ty(F,6). For simplicity, 
in what follows we will associate the notion of maximum 
resistance force specifically with 9 = 0. 

From the condition that the total torque keeps null in 
relation to O, it is not difficult to prove that F = P tan *. 
This relation confirms that * increases with F. It then 
follows from d3j) that / = ||f|| = yP, with \x = tan*, 
where we have introduced /1 to make the link with the 
theory for blocks. Now we look for the maximum value 
of ^ that ensures static equilibrium. Presumably, this 
value depends on the physical properties of the surface, 
but it must not be greater than the geometrical bound 
^> g given by Eq. (Q]), otherwise the point of application of 
the contact force would be out of the contact region. Our 
model concerns situations in which no sliding occurs, so 
that we assume that the rolling will start at some value 
"J for which the surface sinks under the pressure caused 
by the circular body. The perturbation induced by the 



critical force on the surface tends to propagate initially as 
a wave. From this stage on, theories of rolling friction can 
be invoked to explain the motion 0, 0, 0] ; but this is out 
of the scope of this contribution. Given that \I> ^ ty g , the 
maximum resistance force the surface can provide before 
the rolling starts is upper bounded as 



(i) 



/max ^ M -*i 



with 



/' 



tan* g . 



(4) 



This relation [along with Eqs. (HJ and <j2j) ] defines our 
model for the maximum resistance force. Although the 
above relation provides just an upper bound, it furnishes 
a quantitative formula with no free parameter. In addi- 
tion, our result depends only on physical properties that 
are generally known, so that the adequacy of model can 
be experimentally assessed. 



D. Unified view 

There is another remarkable aspect that emerges from 
our simple model, namely, the possibility of unifying 
the results for spheres and cylinders in a single formula. 
Consider the superficial areas S = A-kR 2 (sphere) and 
S = 2-kRL (cylinder). We then define, for future conve- 
nience, the following dimensionless quantity: 



P 



P/S superficial pressure 
n medium constant 



(5) 



The pressure ratio p is clearly a quantity that en- 
capsulates information about both participants of the 
phenomenon, the indenter (the circular body) and the 
medium (the elastic surface). Roughly, it measures the 
local capability of the body in to deform the surface. 

Now, since we are focusing in regimes in which the 
penetration length d is small, one may use Eqs. ([1]) 
and (|4]) to approximate the coefficient of resistance as 
/i = * g = (2d/R) 1/2 . From Eqs. © and flSJ) it is straight- 
forward to show that one may unify the results for the 
sphere and the cylinder into 



^fop 1/(D+1 » 



(0) 



where (3d = (5 — D)^ 2 ~ D ^ 2 . The result for the sphere 
(cylinder) is retrieved with D = 2 (D = 1). Interest- 
ingly, the constant D can be directly associated with the 
symmetry of the contact: While for the sphere the dis- 
tribution of pressure is two-dimensional, for the cylinder 
it is one-dimensional [4j, |7(. It is instructive to imag- 
ine the situation right after the beginning of the contact, 
when the body has just been released on the surface. For 
the sphere, the contour lines of pressure are described by 
concentric circumferences (say, symmetry D = 2) while 
for the cylinder they form concentric thin long stripes 
(symmetry D = 1). This observation further valorizes 
the result ([6]), as it encodes into the resistance coefficient 
information about the symmetry of the contact. 

As far as the parameter (3d is concerned, some com- 
ments are in order. While for the sphere /?2 = 3, for the 



4, in accordance with 

1/2 



cylinder on has that either fi\ 

Ref. i, or / gf ) = [61n(^)] 1/2 (see Ref. (H), the 

latter keeping an explicit dependence on p. Because of 

(2) 

its logarithmic form, j3\ applies as a meaningful physi- 
cal parameter only if p < (L/R) 2 . While this condition 
seems to further restrict the applicability of the model, 
it actually is in full harmony with the premises of our 
approach. To check that, one should note, on one hand, 
that the formula (|2b[) was derived in Ref. [33[ under the 
assumption that L 3> R and, on the other hand, that we 
have previously required that p <C 1, which corresponds 
to the regime of hard clastic media and hence small pene- 

(2) 

tration lengths. It follows, therefore, that /3[ is indeed a 
well behaved parameter in the range p -C (L/R) 2 within 
which our model applies. Furthermore, likewise /?2 and 
P{ , j3j_ does not significantly vary throughout the (in 
principle permissible) several decades for p (see Fig. [3]) , 
so that the unified formulation \i ex p 1 /( D+1 ) reveals to 
be meaningful indeed. 




FIG. 3. Log-log plot for both /i (thick lines) and /3d as a 
function of p — P/nS [see Eq. ©]. The color code is as 
follows: Black for /?2 (sphere), red for j3\ (cylinder, Ref. [J]), 



and blue for f3\ 
L = 5R. 



(2) 



(cylinder, Ref. [33J). In these simulations, 



III. CONCLUDING REMARKS 

The physics of the friction is widely addressed via phe- 
nomenological approaches. The formula fiP is universally 
accepted for both the static and the kinetic friction, with 
/i encapsulating all the unknown physics. The situation 
is not different when the problem of rolling friction is 
concerned. In this paper, we have modeled the horizon- 
tal force a surface applies to a circular body when the 
latter remains at rest while acted upon by an external 
force. To some extent, this resistance force stands for 
the problem of circular bodies as the static friction for 
blocks. To the best of our knowledge, this problem has 
never been addressed so far. 

Our approach is built out from a minimum mechani- 
cal model that conceives the contact of an ideally rigid 



circular body, such as a cylinder or a sphere, with an 
elastic plane surface of known properties. By assuming 
that both the resistance and the normal forces are com- 
ponents of a single contact force C we have derived an 
upper bound for the maximum resistance force that the 
surface can provide before the rolling starts. Our main 
result — composed of Eqs. Q and ^ — shows that the 
maximum resistance force is proportional to the load P 
via a coefficient that encapsulates all the physics of the 
contact. Indeed, /i is shown to depend solely on physical 
properties of the problem, namely, the pressure per unit 
area caused by the circular body, the medium constant, 
defined by a relation between the Young modulus and 
the Poisson ratio, and the symmetry of the contact. In 
particular, the ratio p between the superficial pressure 
(P/S) and the medium constant in units of pressure (k), 
emerges as distinctive figures of merit for the physics of 



Finally, we hope that our approach can motivate fur- 
ther explorations on this basic problem. Naturally, a 
first theoretical improvement would be to study regimes 
described by the JKR and DMT theories [7| , which con- 
sider adhesion and Van der Waals interactions. We be- 
lieve that tests of our model can be carried on with cur- 
rent technologies, in particular within the context of con- 
trolled manipulation of nanoparticlcs [35l . |36| . 
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